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The  new  thermal  rate  equations  were  built  up  by  taking  the  linear  and  non-linear  components  in  the 
temperature  dependences  of  the  Seebeck  coefficient  a,  the  electrical  resistivity  p  and  thermal  conductiv¬ 
ity  k  of  a  thermoelectric  (TE)  material  into  the  thermal  rate  equations  on  the  assumption  that  their  tem¬ 
perature  dependences  are  expressed  by  a  quadratic  function  of  temperature  T.  The  energy  conversion 
efficiency  rj  for  a  single  TE  element  was  formulated  using  the  new  thermal  rate  ones  proposed  here. 
By  applying  it  to  the  high-performance  half-Heusler  compound,  the  non-linear  component  in  the  temper¬ 
ature  dependence  of  a  among  those  of  the  TE  properties  has  the  greatest  effect  on  r\,  so  that  rjho  was 
increased  by  1 1%  under  the  condition  of  T=  510  K  and  AT  =  440  K,  where  rjo  is  a  well-known  conventional 
energy  conversion  efficiency.  It  was  thus  found  that  the  temperature  dependences  of  TE  properties  have  a 
significant  influence  on  rj.  When  one  evaluates  the  accurate  achievement  rate  of  rjexp  obtained  experimen¬ 
tally  for  a  TE  generator,  therefore,  rjexp  should  be  compared  with  rjthe  estimated  from  the  theoretical 
expression  proposed  here,  not  with  rj0,  particularly  when  there  is  a  strong  non-linearity  in  the  tempera¬ 
ture  dependence  of  TE  properties. 
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1.  Introduction 

There  has  been  considerable  interest  during  the  past  10  years  in 
finding  new  materials  and  structures  for  use  in  clear,  highly  effi¬ 
cient  cooling  and  energy  conversion  systems  [1,2].  The  thermo¬ 
electric  (TE)  figure  of  merit  ZT  provides  a  measure  of  the  quality 
of  such  materials  for  applications  and  is  defined  as  ZT  =  a2T/pK, 
where  a  is  the  Seebeck  coefficient,  p  is  the  electrical  resistivity,  k 
is  the  thermal  conductivity  and  T  is  the  absolute  temperature  [1]. 
The  increase  in  ZT  leads  directly  to  the  improvement  in  the  energy 
conversion  efficiency  of  TE  generators  and  in  the  cooling  efficiency 
of  Pelteir  modules  [3].  Much  effort  has  been  made  to  raise  ZT  of  TE 
bulk  materials  for  the  improvement  of  energy  conversion  effi¬ 
ciency  and  there  were  indeed  some  improvements  in  ZT  (for  exam¬ 
ple,  ~ 2  at  700  K  for  the  p-type  TAGS  (GeTe-AgSbTe2)  [4],  about  1.5 
at  700  K  for  the  n-type  (Zr-Hf)NiSn  half-Heusler  compound  [5]) 
and  ~2.2  at  800  K  for  the  n-type  AgPb18SbTe20  [6],  but  even  their 
values  are  not  sufficient  to  improve  dramatically  the  energy  con¬ 
version  efficiency. 

The  energy  conversion  efficiency  rj  has  been  derived  analytically 
using  the  well-known  thermal  rate  equations  at  the  cold  and  hot 
junctions  of  a  TE  element  [3,7].  The  linear  change  in  a  with  temper¬ 
ature  has  been  taken  into  the  conventional  thermal  rate  equations 
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by  Ioffe  [7]  and  the  linear  changes  in  a,  p  and  k  with  temperature 
were  recently  taken  into  them  by  us  [8].  However,  the  non-linear 
components  in  the  temperature  dependences  of  the  TE  properties 
have  been  neglected  in  the  thermal  rate  equations,  although  gener¬ 
ally  they  vary  non-linearly  with  changes  in  temperature  along  a  TE 
element.  The  traditional  expression  for  r\  obtained  from  the  con¬ 
ventional  thermal  rate  equations  should  thus  be  not  correct  ex¬ 
actly.  From  this  reason,  rj  has  been  calculated  numerically  by 
solving  the  differential  equations  governing  TE  properties  by  the 
software  tool  in  which  the  correct  temperature  dependence  of  TE 
properties  is  taken  into  account  [9].  However,  such  a  calculation 
has  never  clarified  analytically  the  degree  to  which  each  tempera¬ 
ture  dependence  of  TE  properties  has  an  effect  on  r\.  To  make  it 
clear,  we  made  an  attempt  to  take  the  difference  between  the  TE 
properties  at  the  hot  and  cold  sides  of  a  TE  element  into  the  ther¬ 
mal  rate  equations  on  the  assumption  that  all  of  the  TE  properties 
are  expressed  by  a  quadratic  function  of  temperature.  The  new 
thermal  rate  equations  were  thus  formulated  by  taking  both  the 
linear  and  non-linear  components  in  the  temperature  dependences 
of  TE  properties  into  the  thermal  rate  equations.  The  reason  that 
the  expression  of  rjlrj0  derived  here  was  applied  to  the  half-Heusler 
compound  [5]  is  that  the  temperature  dependences  of  its  TE  prop¬ 
erties  are  expressed  well  by  a  quadratic  function  of  temperature 
and  its  compound  has  an  excellent  ZT,  where  r\  and  rj0  are  the  en¬ 
ergy  conversion  efficiencies  calculated  using  the  new  and  the  con¬ 
ventional  thermal  equations,  respectively.  r\  was  then  derived 
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Nomenclature 

TE 

thermoelectric 

C2 

coefficient  of  ( Tz  -  T)2  when  k(T)  was  expanded  in  a 

Z 

thermoelectric  figure  of  merit  (K-1) 

power  series  in  (Tz  -  T)  (K  2) 

ZT 

dimensionless  thermoelectric  figure  of  merit 

Qh 

heat  energy  at  the  hot  junction  of  a  TE  element  (W) 

T 

mean  temperature  of  a  TE  element  (given  by 

Qc 

heat  energy  at  the  cold  junction  of  a  TE  element  (W) 

T=(Tf,  +  Tc)/2)  (K) 

Q 

difference  (Qh  -  Qc)  in  the  heat  energy  between  both 

AT 

temperature  difference  between  both  ends  of  a  TE  ele- 

end  junctions  of  a  TE  element  (W) 

ment  (K) 

S 

cross-sectional  area  of  a  TE  element  (m2) 

Tz 

temperature  at  a  position  z  of  a  TE  element  (K) 

L 

length  of  a  thermoelectric  material  of  a  TE  element  (m) 

Th 

temperature  at  the  hot  junction  of  a  TE  element  (given 

Rl 

load  resistance  (Q) 

by  T+AT/2)  (K) 

I 

electric  current  flowing  through  a  TE  element  (A) 

Tc 

temperature  at  the  cold  junction  of  a  TE  element  (given 

by  T-  AT/2)  (K) 

Greek  letters 

I((T) 

thermal  conductance  given  by  I<(T)  =  5k(T)/L  (W/K) 

a 

Seebeck  coefficient  of  a  TE  material  (V/K) 

R(T) 

electric  resistance  given  by  R(T )  =  Lp(T)/5  (Q) 

a  (Tn) 

Seebeck  coefficient  at  the  hot  junction  of  a  TE  element 

(I<(Th)) 

thermal  conductance  averaged  over  a  half  (hot)  side  of  a 

(V/K) 

TE  material  (W/K) 

«(TC) 

Seebeck  coefficient  at  the  cold  junction  of  a  TE  element 

(■ K(TC )) 

thermal  conductance  averaged  over  a  half  (cold)  side  of 

(V/K) 

a  TE  material  (W/K) 

K 

thermal  conductivity  of  a  TE  material  (W/mK) 

( K(T )) 

thermal  conductance  averaged  all  over  a  TE  material 

P 

electrical  resistivity  of  a  TE  material  (Q  m) 

(W/K) 

y 

parameter  defined  as  (1  +ZT)-1/2 

(R(Th)) 

electric  resistance  averaged  over  a  half  (hot)  side  of  a  TE 

P(Th) 

Thomson  coefficient  at  the  hot  junction  of  TE  element 

material  (Q) 

(defined  as  Th(da/dT))  (V/K) 

(R(Tc)) 

electric  resistance  averaged  over  a  half  (cold)  side  of  a 

P(TC) 

Thomson  coefficient  at  the  cold  junction  of  TE  element 

TE  material  (Q) 

(defined  as  Tc(da/dT))  (V/K) 

M 

coefficient  of  (Tz  -  T)  when  a(T)  was  expanded  in  a 

p 

energy  conversion  efficiency  for  ^  0  and  Cj  ^  0  K1 

power  series  in  (Tz  -  T)  (K1) 

and  A2  *  0,  B2  #  0  and  C2  #  0  K  2  (%) 

coefficient  of  ( Tz  -  T)2  when  a (T)  was  expanded  in  a 

ho 

energy  conversion  efficiency  for  =  Cj  =  0  K-1  and 

power  series  in  (Tz  -  T)  (K-2) 

A2  =  B2  =  C2  =  0  K“2  (%) 

Bi 

coefficient  of  (Tz  -  T)  when  p(T)  was  expanded  in  a 

h  exp 

energy  conversion  efficiency  obtained  experimentally 

power  series  in  (Tz  -  T)  (K1) 

for  a  TE  generator  (%) 

b2 

coefficient  of  (Tz  -  T)2  when  p(T)  was  expanded  in  a 

hthe 

energy  conversion  efficiency  estimated  for  a  generator 

power  series  in  (Tz  -  T)  (K-2) 

using  the  theoretical  expression  (%) 

c, 

coefficient  of  (Tz  -  T)  when  k{T)  was  expanded  in  a 
power  series  in  (Tz  -  T)  (K1) 

analytically  on  the  assumption  that  the  terms  dependent  on  the 
temperature  difference  AT  imposed  on  a  TE  element  are  suffi¬ 
ciently  small  compared  with  those  independent  of  AT.  It  is  clarified 
in  this  paper  that  rj/rjo  of  the  high-performance  half-Heusler  com¬ 
pound  [5]  is  affected  significantly  by  the  temperature  dependences 
of  TE  properties.  In  addition,  it  is  investigated  what  effect  each 
temperature  dependence  of  TE  properties  has  on  rjlp0. 

The  purpose  of  this  paper  is  to  investigate  what  effect  the  linear 
and  non-linear  components  in  the  temperature  dependences  of  a, 
p  and  k  of  a  TE  material  have  on  ^  of  a  single  TE  element  and  is 
to  show  to  what  degree  ^  of  a  TE  element  is  affected  by  the  temper¬ 
ature  dependences  of  TE  properties. 

2.  Analysis 

2.1.  Temperature  dependence  of  thermoelectric  properties 

Let  us  consider  a  single  TE  element  in  which  a  TE  parallelepiped 
is  welded  with  two  electrically  conducting  materials,  as  shown  in 
Fig.  1.  When  the  temperature  difference  was  imposed  on  a  TE  ele¬ 
ment,  the  voltage  is  generated  between  both  ends  of  a  TE  element 
and  the  current  flows  through  a  circuit.  The  current  can  then  be 
used  to  power  a  load,  thus  converting  the  thermal  energy  into  elec¬ 
tric  energy.  For  a  majority  of  cases  a  one-dimensional  model  has 
turned  out  to  be  sufficient  to  obtain  quantitatively  correct  perfor¬ 
mance  data  [10].  Assuming  that  the  temperature  difference  be¬ 
tween  the  hot  (Th)  and  cold  (Tc)  junctions  of  a  TE  element  is  AT 
and  the  mean  temperature  is  T,  Th  and  Tc  are  expressed  as  T  +  AT/ 


2  and  T  -  AT/2,  respectively,  at  least  when  the  temperature  within 
a  TE  parallelepiped  changes  linearly  along  its  length  direction. 
Thus,  r\  should  depend  on  the  temperature  dependence  of  the  TE 
properties  of  a  TE  material,  because  generally  the  TE  properties  dif¬ 
fer  appreciably  at  the  hot  and  cold  junctions. 

Generally,  the  Seebeck  coefficient  a(T),  the  electric  resistance 
R(T)  and  the  thermal  conductance  K{T)  of  a  TE  element  vary  with 
changes  in  temperature  T  over  a  working  temperature  range.  The 
difference  in  a  between  both  junctions  has  been  taken  into  the 
conventional  thermal  rate  equations,  but  the  differences  in  R  and 


Fig.  1.  A  single  TE  element  (possessing  the  length  L  and  cross-sectional  area  S) 
sandwiched  between  two  metals  on  which  the  temperature  difference  AT  is 
imposed,  where  Th  and  Tc  are  the  temperatures  at  the  hot  and  cold  junctions,  T  is  the 
mean  element  temperature  and  I  is  a  current  flowing  through  a  circuit  composed  of 
a  TE  element  and  a  load  resistance  RL. 
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I<  between  the  hot  and  cold  sides  have  always  been  neglected.  In 
other  words,  ot  has  been  replaced  with  an  average  value  (a)  over 
a  working  temperature  range,  the  resistances  of  the  hot  and  cold 
sides  have  been  replaced  by  ft/2  so  that  the  Joule  heat  is  divided 
equally  into  the  hot  and  cold  sides  and  I<  has  been  regarded  as  a 
constant  throughout  a  TE  element.  However,  the  linear  changes 
in  oc,  p  and  k  with  temperature  were  recently  taken  into  them  by 
us,  so  that  it  was  clarified  that  the  temperature  dependence  of  p 
has  a  significant  effect  on  27/770  [8].  It  is  thus  examined  here  to  what 
degree  rj  of  a  single  TE  element  is  affected  by  the  linear  and  non¬ 
linear  components  in  the  temperature  dependences  of  TE  proper¬ 
ties.  Here,  let  us  assume  that  a (T),  p(T)  and  tc(T)  of  a  TE  element 
are  expressed  as  a  quadratic  function  of  T  over  a  working  temper¬ 
ature  range.  As  described  in  Appendix,  when  the  temperature 
within  a  TE  parallelepiped  changes  linearly  along  its  length  direc¬ 
tion  (see  Fig.  1 ),  both  (xh  and  ac  at  the  hot  and  cold  junctions  are  ex¬ 
pressed  as 

a„  =  ac  =  a(r)[l-(AT)2A2/4]  (1) 

by  taking  the  Thomson  coefficients  into  the  thermal  rate  equations. 
Subsequently,  the  resistances  ( ft(Th ))  and  (ft(Tc))  averaged  over  a 
half  (hot  or  cold)  side  of  a  TE  parallelepiped  are  expressed  as 

<R(7„))  =  (R(7)/2)[l  +  (AT)E,/4  +  (A7)2B2/12]  (2) 

and 

(R(Tc))  =  (R(7)/2)[l  -  (A7)B,/4  +  (A7)2B2/12]  (3) 

as  a  function  of  AT  (see  Appendix).  The  conventional  theory  has 
been  dealt  as  =  0  K_1  and  B2  =  0  K-2,  but  we  will  treat  here  as 
B  i^O  K1  and  B2  *  0  K“2,  because  (AT^/4  and  (AT)2B2/12  have 
values  of  -0.13  and  0.016  at  AT  =  440  K  in  the  half-Heusler  com¬ 
pound  [5],  whose  values  are  not  negligibly  small  compared  with 
one. 

When  AT  was  applied  to  a  TE  element,  the  thermal  conductance 
(K(T))  averaged  all  over  a  TE  parallelepiped  is  expressed  as 

(R(T)j  =  R(7)[l  -  (AT)2(C2  -  C2)/l 2]  (4) 

as  described  in  Appendix,  where  K(T )  =  Stc(T)/L.  The  reason  that  the 
thermal  conductance  of  a  TE  material  cannot  be  divided  into  two 
parts  of  (K(Th))  and  ( K(TC ))  averaged  over  a  half  side  is  that  the  con¬ 
tinuity  of  heat  conduction  along  a  TE  element  never  permits  us  to 
divide  (K(T))  of  a  TE  material  into  ( I<(Th ))  and  ( I<(TC )).  Even  though 
it  was  divided  into  two  parts,  however,  the  difference  between 
(K{Th))  and  ( I<(TC ))  of  the  half-Heusler  compound  is  small  so  that 
it  is  not  necessary  to  discriminate  between  them,  because  there  is 
little  difference  in  k  between  the  hot  and  cold  junctions.  Therefore, 
it  would  be  allowed  to  express  a  thermal  conductance  of  a  TE  ele¬ 
ment  by  ( K{T ))  alone.  In  the  half-Heusler  compound,  (A TC^)2 1 12 
and  (AT)2C2/12  take  values  of  0.00017  and  0.044  at  AT*  440  K, 
respectively.  The  former  is  much  smaller  than  one,  but  the  latter 
is  not  negligibly  small  compared  with  one.  The  whole  (I<{T))  of  a 
TE  material  depends  on  the  sign  of  C2  but  not  on  the  sign  of  C\. 

As  mentioned  above,  Eqs.  (1  )-(4)  were  derived  by  replacing  dTz\ 
dz  with  A T/L  on  the  assumption  that  a  TE  element  has  a  linear  tem¬ 
perature  profile  along  its  length.  For  this  reason,  these  equations 
may  be  inapplicable  for  the  practical  generators  which  have  a 
strong  non-linearity  in  the  temperature  profile  along  the  length. 
However,  it  is  possible  to  derive  more  accurate  expressions  for 
och,  otc ,  (ft(Tjfj)),  (ft(Tc))  and  (I<(T))  by  substituting  the  practical  tem¬ 


perature  gradient  dTJdz  into  Eqs.  (A.5),  (A.6),  (A.9),  and  (A.13)  in 
Appendix. 

2.2.  r]  derived  by  taking  into  account  the  temperature  dependence  of 
thermoelectric  properties 

Here  we  make  an  attempt  to  take  the  temperature  dependences 
of  TE  properties  into  the  energy  conversion  efficiency  r\.  First  of  all, 
we  derive  the  new  thermal  rate  equations  using  ot(Th),  a(Tc),  (R{Th)), 
(ft(Tc)),  (K(T)),  p(Tft)  and  p(Tc),  where  p(Th)  and  p(Tc)  are  the  Thom¬ 
son  coefficients  defined  as  p {Th)  =  Th(doc/dT)  and  p(Tc)  =  Tc(da/dT) 
[11].  When  AT  is  applied  to  a  TE  element  and  the  Seebeck  and 
Thomson  effects  result  in  the  power  generation,  the  thermal  rate 
equations  at  the  hot  and  cold  junctions  for  a  TE  element  can  thus 
be  expressed  as 


Qh  =  ochThI  -  ( R(7„)>/ 2  +  (K(T))AT 

(5) 

and 

Q.c  =  acTcI  +  (R(7C)>/2  +  (K(T))AT 

(6) 

where  T  is  the  mean  TE  element  temperature  of  (Th  +  Tc)/2,  (I<(T))  is 
the  average  thermal  conductance  of  a  TE  material,  oth  and  ac  are  the 
Seebeck  coefficients  at  the  hot  and  cold  junctions,  and  Qh  and  Qc  are 
thermal  energies  heated  and  cooled  at  the  hot  and  cold  junctions 
[12].  The  first,  second  and  third  terms  in  both  equations  represent 
the  Peltier  and  Thomson  effects,  Joule  heating  and  thermal  conduc¬ 
tion,  respectively.  When  (ft(Th)>  =  (ft(Tc))  =  ft(T)/2  and  (7C(T))  =  JC(T), 
these  equations  correspond  to  the  thermal  equations  reported  by 
Min  et  al.  [11]  By  substituting  Eqs.  (l)-(4)  into  Eqs.  (5)  and  (6), 
the  new  thermal  rate  equations  are  expressed  as 

Qh  =  a(W  ~  (AT)2A2/4]ThI  -  R(T)[  1  +  (AT)B,/4 

+  (A7)2B2/12]/2/2  +  /C(T)AT[1  -  (A 7)2(C2  -  C2)/12]  (7) 

and 

Qc  =  «(T)[  1  -  (A7)2A2/4]7C/  +  R(7)[1  -  (AT)B,/4 

+  (AT)2B2/12]/2/2  +  K'(T)AT[1  —  (A7)2(C2  —  C2)/l  2]  (8) 

at  least  when  there  is  no  heat  transfer  to  or  from  the  surroundings 
other  than  a  TE  element.  After  all,  the  term  of  Ai  is  canceled  out  by 
that  of  the  Thomson  coefficient,  as  already  pointed  out  by  Ioffe  [7], 
so  that  the  terms  of  A2,  Bi,  B2,  Q  and  C2  remain  in  the  thermal  rate 
equations.  If  these  parameters  are  equal  to  zero,  of  course,  these 
thermal  rate  equations  are  reduced  to  the  conventional  thermal  rate 
equations  [7].  The  energy  difference  Q  between  Qh  and  Qc  is  given 
by 

Q  =  Q/i  Qc 

=  a(T)[l  -  (A7)2A2/4]A7J  -  R(7)[l  +  (AT)2B2/12]/2  =  RLI2  (9) 

where  RLI2  is  the  electric  energy  consumed  in  a  load  resistance  RL 
[7].  From  Eq.  (9),  the  current  /  flowing  through  a  circuit  is  obtained 
simply  as 

j  «(7)[1-(A7)2A2/4]A7 
{Rl  +  R(7)[1  +  (AT)2B2/12]} 

The  energy  conversion  efficiency  27  of  a  TE  element  is  defined  as 
QJQh  [7]  and  is  expressed  as 


_ Rf _ 

a(T)[l  -  (A7)2A2/4]7„/  -  R(7)[l  +  (A7)B,/4  +  (A7)2B2/12]/2/2  +  K(7)[1  -  (A7)2(C?  -  C2)/12] 


(11) 
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By  substituting  Eq.  (10)  into  Eq.  (1 1 )  and  setting  RlIR{T)  =  m,  27  is  When  the  terms  of  (AT)2  are  very  small  compared  with  one,  Eq. 

rewritten  as  (13)  is  rewritten  as 


where 

a  =  ZATm[l  -  (AT)2A2/4]2, 
b  =  [1  -  (A T)2(C2  -  C2)/12][m  +  1  +  (AT)2B2/U)}2 

and 

c  =  Z[  1  -  (AT)2A2/4]2[Tfl(m  +  1  -  (AT)2B2/12) 

-  (AT/2)(1  +  (A7)B,/4  +  (AT)2B2/12)] 

where  Z  =  y.{T)2 1 R(T)K(T).  By  setting  dr//dm  =  0,  the  quadratic  equa¬ 
tion  of  m  for  the  maximum  /;  is  obtained  as 

(1  -m2)[l  -(AT)2(C2-C2)/12]+ZT-Z(AT)2B1/8 

+  (1  +  ZT/2)(AT)2B2/6  -  Z(AT)2A2/2  =  0  (13) 


200  400  600  800 

T(  K) 

Fig.  2.  Temperature  dependences  of  (a)  a,  (b)  p  and  (c)  k  for  the  n-type  half-Heusler 
compound. 


m2  =  r2{l  +  [(1  -  y2)(A7-)2/12][6A2  -  (3B,/27) 

+  (1  +  y2)B2/(l  -  y2)  +  (C2  -  C2)]}  (14) 

where 


y  =  (l+ZT)-1/2  (15) 

When  the  second  term  in  the  brace  in  Eq.  (14)  is  very  small 
compared  with  one,  m  is  expressed  approximately  as 

m  =  r’{l  +  (1  -  y2)(A7)2/24][6A2  -  (3B,/2T) 

+  (1  +  y2)B2/(l  -  y2)  +  (C2  -  C2)]}  (16) 


By  substituting  Eq.  (16)  into  Eq.  (12),  moreover,  27  is  expressed 
by  a  similar  approximation  as 


n  =  n  oi  i  + 


(i  +  y)y(A7)2r 
2(yTc  +  Th) 


-a2  + 


(i-y)Bi  ,  (-b2  +  c2-c2) 

4(1  +y)T 


+ 


(17) 


where  rj0  is  the  original  energy  conversion  efficiency  derived  from 
the  conventional  thermal  rate  equations  and  is  expressed  as  [7] 


(l-y)AT 
/0  (yTc  +  Th) 


(18) 


When  a,  R  and  K  are  expressed  by  a  quadratic  function  of  T,  the 
term  of  A\  disappears  but  the  terms  of  other  parameters  remain  in 
27.  Of  course,  r\  is  reduced  to  p0  for  =  Ci  =  0  K-1  and 
A2  =  B2  =  C2  =  0 1<-2.  It  is  thus  found  that  27/270  depends  strongly  on 
the  sign  and  magnitude  of  parameters  of  A2 ,  B2 ,  CA  and  C2.  Eq. 
(17)  is  found  to  provide  an  information  enough  to  investigate  the 
effect  of  the  temperature  dependences  of  a,  p  and  k  on  27,  although 
it  is  derived  approximately  from  the  thermal  rate  equations. 

However,  it  should  be  noted  here  that  even  when  the  present 
approximation  method  is  unavailable,  it  is  possible  to  estimate 
the  accurate  27  numerically  without  approximation.  That  is,  27  can 
be  estimated  exactly  by  substituting  the  value  of  m  obtained 
numerically  from  Eq.  (13)  using  the  values  of  various  parameters 
into  Eq.  (12). 


3.  Results  and  discussion 

3.1.  Parameters  ofAh  A2,  Bh  B2,  Cj  and  C2for  the  high-performance 
half-Heusler  compound 

Figs.  2  and  3  show  the  temperature  dependences  of  a,  p ,  k  and  Z 
for  the  high-performance  n-type  (Zr0.5Hfo.5)Tio.5NiSno.998Sb0.oo2 
compound  [5],  which  has  an  extremely  high  value  of  ZT  =  1.50  at 
700  K.  The  TE  properties  of  this  compound  can  be  expressed  well 
by  a  quadratic  function  of  (T-  510),  as  shown  in  Fig.  2.  Tz  used 
in  Appendix  will  be  replaced  by  T  hereafter,  because  T  used  until 
this  time  is  replaced  with  a  constant  value  of  T=  510  K.  The  param¬ 
eters  of  Ai,  A2,  B^,  B2i  Ci  and  C2  which  are  coefficients  of  (T-  510) 
and  (T-  510)2  were  determined  by  the  least  square  fitting  to  the 
experimental  values  of  a,  p  and  k  for  the  half-Heusler  compound 
in  the  temperature  range  from  290  to  730  K,  as  listed  in  Table  1. 
In  this  compound,  Ai,  B2,  Ci  and  C2  are  positive  in  sign  but  A2 
and  B 1  are  negative  in  sign.  Of  course,  the  values  of  these  parame¬ 
ters  would  change  significantly  with  the  chemical  composition  of 
half-Heusler  compounds  [5,13-15].  The  Z  value  calculated  using 
these  parameters  reproduces  well  their  experimental  values,  as 
shown  in  Fig.  3a.  This  good  agreement  between  the  experimental 
values  and  calculated  curve  indicates  that  the  present  approxima¬ 
tion  is  valid  for  this  half-Heusler  compound  employed  here.  The 
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Fig.  3.  Temperature  dependences  of  Z  and  ZT  measured  for  the  n-type  half-Heusler 
compound,  where  the  solid  curve  of  Z  was  drawn  using  the  expressions  obtained  by 
the  least  square  fitting  to  the  temperature  dependences  of  a,  p  and  k  and  that  of  ZT 
was  drawn  to  pass  through  the  experimental  values. 


Table  1 

The  parameters  of  A,  B  and  C  estimated  at  T=510K  from  the  temperature 
dependences  of  the  thermoelectric  properties  for  the  half-Heusler  compound  (Ref. 
[5]). 


temperature  dependences  of  TE  properties  are  expressed  well  by  a 
quadratic  function  of  temperature  in  some  high-performance  TE 
materials;  for  example,  AgPbi8SbTe2o  [8]  and  Bi-Te  [16]. 

The  Thomson  coefficients  p(Th)  and  p(Tc)  for  the  half-Heusler 
compound  were  calculated  using  the  relations  p(Th)  =  Th(da/dT)T=Th 
and  p(Tc)  =  Tc(da/dT)T=Tc.  As  a  result,  p (Th)  at  Th  =  730  K  and  p(Tc)  at 
TC  =  290K  were  138.1  and  -129.8  pV/K,  respectively.  Their  abso¬ 
lute  values  are  close  to  p  =  129pV/K  estimated  theoretically  by 
Ioffe  [7]  using  a  non-degenerated  model.  It  indicates  that  the  va¬ 
lence  band  of  the  half-Heusler  compound  may  indeed  form  a 
non-degenerate  band. 

3.2.  Dependence  of  rj/rj0  on  AT,  the  values  of  parameters  and  ZT 

The  dependence  of  27/270  on  AT  was  calculated  as  a  function  of 
AT  by  substituting  the  values  of  parameters  (listed  in  Table  1)  for 
the  n-type  half-Heusler  compound  into  Eq.  (17),  as  shown  in 
Fig.  4.  Consequently,  27/270  increases  quadratically  with  an  increase 
of  AT  and  reached  a  great  value  of  1.108  at  AT  =  440  K.  Then,  the 
degrees  of  contribution  from  the  first,  second  and  third  terms  in 
the  square  bracket  of  Eq.  (17)  to  27/270  are  0.154,  -0.006  and 


Fig.  4.  (a)  ri/rjo  calculated  as  a  function  ofX(=f?!  and  Ci)  and  (b)  rj/r]0  as  a  function  of 
X(=A2,  B2  and  C2)  under  the  temperature  condition  of  T  =  510  K  and  AT  =  440  K  for 
the  n-type  half-Heusler  compound. 


Fig.  5.  rjlrjo  calculated  as  a  function  of  AT  at  T=  510  K  for  the  n-type  half-Heusler 
compound.  The  arrows  in  the  figure  denote  the  value  of  each  parameter. 


-0.040,  respectively.  It  is  thus  found  that  the  degree  of  increase 
in  27/270  depends  strongly  on  the  sign  and  magnitude  of  A2.  If  the 
values  of  these  parameters  remain  unvaried  with  changes  in  AT, 
therefore,  27/270  would  be  enhanced  significantly  with  further  in¬ 
crease  in  AT. 

Next,  the  dependence  of  27/270  on  X(=A2,  Bi,  B2,  Ci  and  C2)  was 
calculated  at  T=  510 1<  as  a  function  of  X  using  Eq.  (17)  under  a  con¬ 
dition  of  AT  =  440  K  for  the  half-Heusler  compound  [5]  and  drawn 
as  a  function  of  X  in  Fig.  5.  As  evident  from  Eq.  (17),  27/270  changes 
linearly  with  an  increase  of  A2,  Bi,  B2  and  C2  but  increases  quadrat¬ 
ically  with  an  increase  of  Ci,  so  that  the  dependences  of  27/270  on  the 
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parameters  change  entirely  according  to  their  sign  and  magnitude. 
It  is  seen  from  the  figure  that  /7/r/o  varies  significantly  with  changes 
in  A2,  particularly  at  large  AT,  as  mentioned  above.  The  sign  and 
magnitude  of  A2  have  a  much  more  significant  effect  on  rjlrj0  than 
those  of  other  parameters.  The  error  due  to  the  approximation  for 
the  increment  ( rj/rio  -  1)  in  r]lrj0  is  about  6%  at  AT  =  440  K,  so  that 
Eq.  (17)  is  found  to  hold  well  even  for  a  great  temperature  differ¬ 
ence  of  AT  =  440  K.  Therefore,  the  present  results  would  be  almost 
unvaried  basically  even  if  rj/rjo  was  derived  numerically  or  exactly 
from  Eqs.  (12)  and  (13). 

Subsequently,  r\  and  rj0  were  calculated  at  T=  510  K  as  a  func¬ 
tion  of  ZT  for  AT  =  440  K  and  their  curves  were  drawn  as  a  function 
of  ZT  in  the  range  from  ZT  =  0  to  3  in  Fig.  6a,  where  r\  was  calculated 
from  Eq.  (17)  using  the  parameters  listed  in  Table  1.  Here  the 
change  in  ZT  means  a  change  in  Z.  Both  rj  and  r]0  tend  to  increase 
monotonically  with  an  increase  of  ZT.  r\  is  somewhat  higher  than 
r\ 0  all  over  the  range.  Under  the  condition  of  T=510K  and 
AT  =  440  K,  however,  rj  is  approximately  11%  at  ZT  =  0.945  higher 
than  rjo  in  the  half-Heusler  compound.  However,  rj/rj 0  tends  to  de¬ 
crease  slightly  with  an  increase  of  ZT,  as  shown  in  Fig.  6b.  In  any 
case,  the  fact  that  r\  is  larger  than  rj0  indicates  that  when  one  eval¬ 
uates  the  accurate  achievement  rate  of  rjexp  obtained  experimen¬ 
tally  for  a  TE  generator,  therefore,  rjexp  should  be  compared  with 
rjthe  estimated  theoretically  from  Eq.  (17),  not  with  the  conven¬ 
tional  r\  0. 

The  estimation  of  rjlr]0  using  the  temperature  dependences  of 
the  TE  properties  evaluated  at  the  average  temperature  T  would 
thus  provide  the  exact  results,  at  least  as  long  as  the  temperature 
dependence  of  TE  properties  was  expressed  as  a  quadratic  function 
of  temperature  and  the  real  temperature  gradient  dTJdz  was  em¬ 
ployed.  Here  the  real  dTz\dz  corresponds  to  the  temperature  gradi¬ 
ent  of  the  temperature  distribution  obtained  by  the  least  square 
fitting  of  a  polynomial  to  the  practical  temperature  profile  of  a 
TE  element. 


Fig.  6.  (a)  r\  and  (b)  rjlrjo  calculated  at  T=510K  as  a  function  of  ZT  under  a 
temperature  difference  of  AT  =  440  K  for  the  n-type  half-Heusler  compound. 


3.3.  ZeffT  obtained  by  ascribing  the  increase  in  r\/r}0  to  the  increase  in 
the  intrinsic  ZT 


The  effective  ZeffT  for  a  TE  element  was  calculated  by  ascribing 
the  increase  in  rilf]0  due  to  the  parameters  of  A2 ,  Bi,  B2 ,  Ci  and  C2 
to  the  increase  in  the  intrinsic  ZT.  For  simplicity,  Eq.  (17)  is  rewrit¬ 
ten  as 

T=1+Y  C9) 

do 

where 


y(l  +  y)(AT)2T 
2  (yTc  +  Th) 


—A2  + 


(i-r)fli 

4(1+7) 


(-B2  +  C?-C2)' 

6 


(20) 


where  y  =  (1  +  ZT)  1/2.  If  we  attribute  this  increase  in  rj/rjo  to  the  in¬ 
crease  in  the  intrinsic  ZT,  r\  is  expressed  as 


„  q-TW)Ar 

1  ( yeffTc  +  Th ) 


(21) 


using  the  effective  y,,ff  =  ( 1  +  ZeffT )  1/2.  Substituting  Eqs.  (18)  and 
(21)  into  Eq.  (19),  is  expressed  as 

(i  —  y)(yTc  +  Th)Y 


yeff  =  y- 


2  T 


(22) 


Substituting  Eq.  (20)  into  Eq.  (22),  the  effective  ZeJ jT is  expressed 
approximately  as 


ZeffT  =  ZT{  1 


(AT)2 


—A? 


(i  -  y)B\  (C2-b2-c2) 

4(1  +  y)T  +  6 


(23) 


Fig.  7.  (a)  The  effective  (Ze//T)  calculated  at  T=  510  K  as  a  function  of  the  intrinsic  ZT 
using  Eq.  (23)  by  ascribing  the  increase  in  rj  due  to  the  parameters  to  the  increase  in 
the  intrinsic  ZT  and  (b)  (ZeJjT)/(Zr)  calculated  as  a  function  of  the  intrinsic  ZT  using 
Eq.  (23). 
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using  the  intrinsic  ZT,  A2,  B2 ,  Q  and  C2.  This  expression  shows 
that  ZeffT  becomes  larger  or  smaller  than  the  intrinsic  ZT  according 
to  the  sign  and  magnitude  of  A2,  Bi,  B 2,  Ci  and  C2.  In  order  to  inves¬ 
tigate  the  effect  of  the  increase  in  27/270  due  to  these  parameters  on 
ZeffT,  therefore,  ZeffT  and  ZeffT I  ZT  were  calculated  as  a  function  of  ZT 
at  T=  510  K  by  substituting  these  parameters  (listed  in  Table  1)  into 
Eq.  (23)  under  the  conditions  of  AT  =  240,  340  and  440  K.  Of  course, 
ZeffT  is  well-defined  by  Eq.  (23)  even  at  AT=440K,  because  the 
approximation  used  in  deriving  Eq.  (17)  holds  well  even  for  a  large 
value  of  AT  =  440  K,  as  mentioned  earlier.  As  shown  in  Fig.  7a,  ZeffT 
increases  almost  linearly  with  an  increase  of  ZT  and  the  gradient 
tends  to  increase  with  an  increase  of  AT.  The  ratio  of  (ZeJyT)  to 
(ZT)  is  enhanced  significantly  at  larger  AT,  but  it  tends  to  decrease 
slightly  and  almost  linearly  with  an  increase  of  ZT,  as  shown  in 
Fig.  7b.  The  ratio  reached  a  great  value  of  1.16  at  ZT  =  0.945  under 
the  condition  of  T=  510  K  and  AT  =  440  K.  It  indicates  that  ZeffT  of 
a  TE  element  is  enhanced  significantly  even  by  the  temperature 
dependences  of  TE  properties.  It  was  thus  found  that  the  tempera¬ 
ture  dependences  of  TE  properties  have  a  significant  influence  on 
the  energy  conversion  efficiency  27  of  a  TE  element,  particularly  at 
large  AT. 


When  the  Seebeck  coefficient  a  changes  linearly  with  T,  a (Tz)  at 
Tz  is  then  expressed  as 

a(Tz)  =  a(T)[  1  +  (Tz  -  T)A ,  +  (Tz  -  T)2A2 ]  (A.2) 

where  A,  =  (da(rz)/dTz)/a(I)  and  A2  =  (d2a(Tz)/df) /2a(T).  As 
shown  in  Subsection  2.1,  ah  and  otc  at  the  hot  and  cold  junctions 
are  expressed  as 

%h  =  oc(Th)  -  p(Th)AT/(2Th)  (A.3) 

and 

uc  =  oc(Tc)  +  /?(TC)AT/  (2TC)  (A.4) 

where  p(Th)  and  p(Tc)  are  the  Thomson  coefficients  defined  by  the 
following  expressions: 

P(Th)  =  Thldcc(Tz)/dTz]Tz=Th  =  Th[(da(Tz)/dz)/(dTz/z)]z=L 

=  Th(A,+ATA2)a(T)  (A.5) 

and 

P(Tc)  =  T  c[doc(T z)  /  dT z\Tz=Tc  =  Tc[(da(Tz)/dz)/(dTz/z)]z=0 

=  TC(A !  -ATA2)oc(T)  (A.6) 


4.  Conclusions 

The  new  thermal  rate  equations  were  built  up  by  introducing 
both  the  linear  and  non-linear  components  in  the  temperature 
dependences  of  a,  p  and  2c  of  a  TE  material  on  the  assumption  that 
they  are  expressed  as  a  quadratic  function  of  temperature.  The  dif¬ 
ferences  in  a  and  p  between  the  hot  and  cold  junctions  of  a  TE  ele¬ 
ment  were  then  taken  into  the  thermal  rate  equations.  The  term  of 
Ai  (coefficient  of  the  linear  temperature  dependence  of  a)  disap¬ 
pears  but  the  terms  of  other  parameters  (A2,  B i,  B2,  Cj  and  C2)  re¬ 
main  in  the  thermal  rate  equations.  Of  course,  when 
£>i  =  Ci  =  0  K-1  and  A2  =  B2  =  C2  =  0  K-2,  the  new  thermal  rate  equa¬ 
tions  are  reduced  to  the  conventional  thermal  rate  ones.  27/270  for  a 
single  TE  element  were  calculated  as  functions  of  AT,  X  and  ZT 
using  the  new  thermal  rate  equations,  where  X  =  A2,  Bi,  B2,  Cj 
and  C2.  Among  the  temperature  dependences  of  TE  properties, 
the  non-linear  component  (A2)  in  the  temperature  dependence  of 
a  has  the  strongest  effect  on  27/270,  so  that  the  degree  of  increase 
in  27/270  of  the  half-Heusler  compound  reached  11%  under  the  con¬ 
dition  of  T=  510  K  and  AT  =  440  K.  When  one  ascribed  this  increase 
in  27/270  due  to  A2,  B 1,  B2,  Ci  and  C2  to  the  increase  in  the  intrinsic  ZT, 
the  effective  ZeffT  of  the  half-Heusler  compound  was  found  to  be 
approximately  16%  at  ZT  =  0.945  higher  than  the  intrinsic  ZT  under 
the  condition  of  T=  510  K  and  AT  =  440  K.  It  was  clarified  that  the 
energy  conversion  efficiency  depends  strongly  on  the  temperature 
dependences  of  a,  p  and  2c. 

When  one  evaluates  the  accurate  achievement  rate  of  rjexp  ob¬ 
tained  experimentally  for  a  TE  generator,  therefore,  rjexp  should 
be  compared  with  rjthe  estimated  theoretically  from  Eq.  (17),  not 
with  the  conventional  rj0. 

Appendix  A 

Let  us  assume  that  the  Seebeck  coefficient  a,  the  electric  resis¬ 
tance  R  and  the  thermal  conduction  I<  of  a  TE  parallelepiped  in¬ 
crease  or  decrease  linearly  with  an  increase  of  T.  Here  we  will 
show  that  (R  (Th))  or  (R  (Tc))  averaged  over  a  half  (hot  or  cold)  side 
of  a  TE  parallelepiped  with  length  L  and  the  resultant  (I<(T))  aver¬ 
aged  all  over  a  TE  parallelepiped  are  expressed  as  a  function  of  AT. 
When  we  assume  that  the  temperature  in  a  TE  parallelepiped 
changes  linearly  with  z  along  its  length  direction,  as  shown  in 
Fig.  1,  the  temperature  Tz  at  a  position  z  is  expressed  as 

TZ  =  T  +  (z-  L/2)AT /L  (A.l) 


Substituting  Eq.  (A.2)  into  Eqs.  (A.3)-(A.6),  the  term  of  Ai  disap¬ 
pears  in  their  equations  but  that  of  A2  remains  so  that  both  och  and 
ac  are  rewritten  as 

ah  =  occ  =  a(T)[l  -  (AT)2A2/4]  (A.7) 

It  is  thus  found  that  the  relation  oth  =  occ  holds  even  for 
A2¥^  0 1<-2  in  the  present  approximation. 

Subsequently,  when  the  electrical  resistivity  p  changes  linearly 
with  T,  p(Tz)  at  Tz  is  then  expressed  as 

p(Tz)  =  p(T)[  1  +  (Tz  -  T)B ,  +  (Tz  -  T)2B2 ]  (A.8) 

Therefore,  the  average  (R  (Th))  of  the  upper  (hot)  side  (corre¬ 
sponding  to  a  half  of  a  TE  parallelepiped)  is  given  by 

mT‘))-t^dz=r^w,dT:  (a-9> 

Using  the  relation  dTz\dz  =  A T/L,  (R(Th))  is  obtained  as 

(R(Th))  =  (R(T)/2)(1  +  (A7)fii/4  +  (A7)2B2/12)  (A.10) 

In  the  same  way,  the  average  (R(TC))  is  also  derived  as 

(R(TC))  =  (R{T)/ 2)(1  -  (AT)B,/4  +  (AT)2B2/12)  (A.ll) 

Moreover,  if  the  thermal  conductivity  k  changes  linearly  with 
temperature,  k(Tz)  at  Tz  is  then  expressed  as 

k(Tz)  =  k(T)[  1  +  (Tz  -  T)C ,  +  (Tz  -  T)2C2]  (A.  12) 

When  there  is  a  temperature  difference  AT  along  a  TE  element, 
therefore,  the  resultant  (I<(T))  of  a  TE  parallelepiped  is  given  by 


1 

Wn) 


=  i  [L-^-dz=\  f 

$  Jo  k(T)z  S  JTc 

—  f 

(T)ATJ. 


/>AT/2 


1  fTh  1  dz 
K(Tz)df 

dx 


dT, 


Sk(T)AT y_|4r/2  (1  +xC\  +x2C2) 


(A.13) 


=  2C2SKL(T)ATAiM[A  +Ar/(Ci/C2  -2Ai)]/[l  - AT/(Ci/C2 -2Ai)]| 
— In|[l  +  A T/(C,  /C2  +  2Ai)]/[l  -  AT/C, /C2  +  2Ai)]|} 


(A- 14) 


where  A  =  (1  -C2/4C2)1/2  and  i  denotes  the  imaginary  number. 
When  |AT/(Ci/C2  +  2Ai)|  «  1  and  |AT/(Ci/C2  -  2Az)|  «  1,  l/(/C(T)>  is 
expanded  in  a  power  series  in  AT  and  the  resultant  (I<(T))  is 
expressed  as 
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1975 


(K(T))  =  K(T) [1  -  (A T)2(C]  -  C2)/12]  (A.  15) 

where  JC(T)  =  Sk(T)/L 
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